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\^ " Abstract 

Gv . 

■ Hamiltonians of a wide-spread class of Gj n „-invariant nonlinear quantum models, in- 

t-h . eluding multiboson and frequency conversion ones, are expressed as non-linear functions 

of si (2) generators. It enables us to use standard variational schemes, based on si (2) gen- 
eralized coherent states as trial functions, for solving both spectral and evolution tasks. 
"^n ■ In such a manner a new analytical expression is found for energy spectra in a mean-field 

CN ! approximation which is beyond quasi-equidistant ones obtained earlier. 

i PACS numbers: 03.70; 02.20 

> : 
\o . 

1 Introduction 

o . 

For last decades a great attention has been paid to developments of both exact and approximate 
techniques to solve and examine different dynamical problems for quantum strongly coupled 
Qh! systems whose interaction Hamiltonians are expressed by nonlinear functions of operators de- 

scribing subsystems (see, e.g., [1-9] and references therein). However, as a rule, such techniques 
either are adapted for treating special forms of model Hamiltonians and initial quantum states 
[1-5,7-9] or require lengthy and tedious calculations (as it is the case, e.g., for the algebraic 
Bethe ansatz [6]). 

Recently, a new universal Lie-algebraic approach has been developed [10-13] to get exact so- 
lutions of both spectral and evolution problems for some nonlinear quantum models of strongly 
coupled subsystems having symmetry groups Gi nv . It was based on exploiting a formalism of 
polynomial Lie algebras g p d as dynamic symmetry algebras g DS of models under study, and, 
besides, generators of these algebras g p d can be interpreted as Gi nv - invariant "essential" collec- 
tive dynamic variables in whose terms model dynamics are described completely. Specifically, 
this approach enabled us to develop some efficient techniques for solving physical tasks in the 
case of g DS = sl p d{2), when model Hamiltonians H are expressed as follows 

H = aV Q + gV+ + g*V- + C, [V a , C] = 0, V. = (V+) + , (1.1) 

where C is a function of model integrals of motion Ri and Vo, V± are the sl p d(2) generators 
satisfying the commutation relations 

[v , \4] = ±y±, [v., v+] = Mv + 1) - MVo), 

n 

MV ) = A 1 [[(V + X l ({R J })) (1.2) 

i=l 
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The structure polynomials ip n (Vo) depend additionally on {Ri,i = 1,...}, and their exact 
expressions for some wide-spread classes of concrete models were given in [10-12]. 

All techniques [10-13] are based on using expansions of most important physical quantities 
(evolution operators, generalized coherent states (GCS), eigenfunctions etc.) by power series 
in the sl p d(2) shift generators V± and on decompositions 

L(H) = J2 L (ih]), (V + V_-i; n (V ) = -MRo)\w]) = (1-3) 

® 

of Hilbert spaces L(H) of quantum model states in direct sums of the subspaces £([/i]) which 
are irreducible with respect to joint actions of algebras sl pd {2) and symmetry groups Gi nv 
and describe specific "s/ p d(2)-domains" evolving independently in time under action of the 
Hamiltonians (1.1); [lo] are lowest weights of L{[lj\) : ip n (lo) = and other quantum numbers 
— 1, . . . are eigennumbers of operators Ri. Then, using restrictions of Eqs. (1.1)-(1.2) on 
L([/j]), one can develop simple algebraic calculation schemes for finding evolution operators 

oo f-1 

U H (t)= E V{u(v ;t), V+ f = v!([MV )] {f Y\ [M^] if) = UM^-r), (1.4a) 

/=-oo r=0 

amplitudes Q v (Ef) of expansions 

\E f ) = A f H(V + - 4)\[k\) =Y,Qv(Ef)\M,v) (1.46) 

3 v 

of energy eigenstates \Ef) in orthonormalized bases : Vo|[Zj];i>) = (lo + v )\[h}', v)}) and 

appropriate energy spectra {£"/} of bound states [11,13]. (In fact, the factorized form of \Ef) 
given by the first equality in (1.4b) realizes an efficient modification of the algebraic Bethe ansatz 
[6] in terms of collective dynamic variables related to the sl p d(2) algebras [11,13].) In the paper 
[12] some explicit integral expressions were found for amplitudes Q V (E), eigenenergies {E a } 
and "coefficients" u(vo]t) of evolution operators Unit) with the help of a specific "dressing" 
(mapping) of solutions of some auxiliary exactly solvable tasks with the dynamic algebra sl{2). 

However, all exact results obtained do not yield simple working formulas for analysis of 
models (1.1) and revealing different physical effects (e.g., a structure of collapses and revivals 
of the Rabi oscillations [2,8], bifurcations of solutions [5] etc.) at arbitrary initial quantum states 
of models. Therefore, it is necessary to develop some simple techniques, in particular, to get 
some closed, perhaps, approximate expressions for evolution operators, energy eigenvalues and 
wave eigenfunctions, which would describe main important physical features of model dynamics 
with a good accuracy (cf. [5,8,9]). Below we examine some possibilities along these lines for 
models (1.1)-(1.2) by means of reformulating them in terms of the formalism of the usual sl(2) 
algebra and developing variational schemes corresponding to quasiclassical approximations for 
original models by analogy with developments [5,14-16]. 

2 A reduction of linear sl p d(2) problems to non-linear sl(2) 
ones 

We can reformulate models (1.1)-(1.2) in terms of s/(2) generators using an isomorphism of 
the sl pd (2) algebras to extended enveloping algebras U,p(sl(2)) of the familiar algebra sl(2). 
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This isomorphism is established via a generalized Holstein-Primakoff mapping given on each 
subspace £([/«]) as follows [10,11] 

y = v - i oT j, y + = ^„_ 2 (y )]-V2, 0n _ 2 (y o) = ^y^ViXyoV Y ~ = 

[Y , Y ± ] = ±Y ± , [Y_, Y + ] = T 2Y (2.1) 

where Y a are the si (2) generators, =Fj are lowest weights of si (2) irreducible representations 
realized on subspaces £([/;]) and ip 2 (x) = (j + l — x) are quadratic structure functions 

ijj n (x) = ip 2 (x) of sl(2) (hereafter upper/lower signs corresponding to the su(2)/su(l, 1) algebras 
are chosen for finite/infinite dimensions d([li]) of the spaces L([/j])). 

Note that, by definition, functions (/> n -2(Yo) on spaces Z/([Z^]) can be chosen as polynomials 
of (n — 2)-th degree in Y . For example, substituting 

^j 3 (V ) = ^(2V + R 2 - R 1 )(2V + R 1 + R 2 )(-V + R 2 + 1), 

k = h = k,l 2 = IM + ^£ k = 0, ±1, ±2, s = 2j = d([h}) - 1 = 0, 1, 2, ... (2.2) 

for three-boson models [11-13] 

H tb = Ld 1 a'la 1 + uo 2 a\a 2 + ^303 a 3 + #(a+a^)a 3 + g*(a 1 a 2 )a 3 , (2.3a) 

V = (iVi + N 2 - N 3 )/3, V + = (aj"a£)a 3 , a = u 1 +u) 2 -u) 3 , N i = ala i) 

2C = Riioj! - u 2 ) + R^ +lo 2 + 2u 3 ), Ri = N ± - N 2 , 3R 2 = N, + N 2 + 2N 3 (2.36) 

we get 4>i(Yo) — Y + j + \k\ + 1. Similar expressions can be found for 4>i(Yq) in the cases of 
the point-like Dicke and the second harmonic generation models taking appropriate expressions 
from [12]. 

Then, restrictions of Hamiltonians (1.1) on £([/«]) may be re- written in terms of Y a as 
follows 

H u = aY + Y + ~g(Y ) + g + (Y )Y_ + C, ~g(Y ) = g^ n _ 2 (Y ), C = C + a(±j + Z„), (2.4) 

Evidently, this form corresponds to generalizations of semi-classical (linear in sl(2) genera- 
tors) versions of matter-radiation interaction models [8,9,12] by introducing operator (intensity- 
dependent) coupling coefficients g(Y ) (cf. [3,7]). Emphasize, however, a collective (not associ- 
ated with a single subsystem) nature of operators Y a in Eq. (2.4) (cf. [9]); therefore, dynamic 
variables Y a correspond to a non-standard quasiclassical approximation (when g(Y ) = const 
in Eq. (2.4)) of original models as it follows, e.g., from a direct comparison of such an ap- 
proximation with standard (when creation/destruction operators of one mode are replaced by 
c-numbers) semiclassical limits for the model (2.3). 

If n = 2, then n _ 2 (Yo) = 1, sl p d{2) = sl{2), l = ±j, and the formalism of GCS related to 
the SL{2) group displacement operators 

S Y (£ = re ie ) = exp(£r+ - £*Y1) = exp[t(r)e ie Y + ] exp[-2 In c(r)Y ] exp[-t(r)e _ *y_], (2.5) 
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(t(r) = tanr/ tanh r,c(r) = cosr/ coshr for su(2)/ su(l, 1)) yields a powerful tool for solving 
both spectral and evolution tasks [16]. 

Specifically, in this case, using the well-known sl(2) transformation properties of operators 
Y a under the action of SV(£) [16]: 

S Y (t)Y + Sy(tf = Y + (0 = [c(r)] 2 Y + ± e^\s(2r)Y, - e -*[ a (r)] 2 Yl], Yl(fl = (y + (0) f , 

5y(Olo^(O t = >o(0 = c(2r)F - ^[e*y+ + e~ ie Y_, s(r) = stnr/stnhr, (2.6) 
Hamiltonians Hy.] can be transformed into the form 

H U (0 = Sy(0%]Sy(O f = + F oA)(a, <?! + Y + A + (a, g; £) + F_A;(a, g; £) (2.7a) 

At the values £o — 2§| arc ^ an ¥ ^ or SM (2) and £o — 2^| arc ^ an ^¥ ^ or SM (1> 1) °f the parameter 
£ one gets A + (a, = 0, and the Hamiltonian H^(^) takes the form 

# [y (£o) = c + y V« 2 ±%l 2 ( 2 - 7& ) 

which is diagonal on eigenfunctions |[/i];u) = ^(j, i>)(y+) u |[Zj]; v = 0), N~ 2 (j,v) = v\(2j)\/(2j — 
v)\ for sw(2) and N~ 2 (j,v) = vW(2j + v)/T(2j) for sw(l, 1). Therefore, original Hamiltonians 
if[y have the eigenenergies 

^(M; Co) = C + {T3 + ^)V« 2 ±4b| 2 (2.8a) 

and eigenfunctions 

\U,v^o)=Sy^ ) t \[h};v;) (2.86) 

Similarly, when sl p d{2) = si (2), operators Sy(£(i)) are "principal" parts in the evolution 
operators Unit) = exp(i(j)(t)Y )SY(^(t)) with c-number functions 4>(t),£,(t) being determined 
from a set of non-linear differential equations corresponding to classical motions [16,17]. 

However, for arbitrary degrees n of ^ n (Vo) Hamiltonians (2.4) are essentially non-linear in 
si (2) generators Y a , and, therefore, the situation is very changed. Particularly, in general cases 
it is unlikely to diagonalize if[y with the help of operators Sy(£) since analogs of Eq. (2.7a) 
on multi-dimensional spaces £([/»]) 

H [h] (0 = Sy(0H [k] S Y (0 j = aY (t) + Y + (0~g(Yo(0) + ~9 + (Y (0)Y- + C (2.9) 

contain (after expanding them in power series) many terms with higher powers of Y± [13]. 

Nevertheless, the formalism of the SL(2) group GCS |[/j];t>;£) = Sy(C,y\[li\;v;) [16] can 
be an efficient tool for analyzing non-linear models [5,11,14-16], in particular, for getting ap- 
proximate analytical solutions. Specifically, a simplest example of such approximations was 
obtained in [11] by mapping (with the help of the change V a — > Y a ) Hamiltonians (1.1) by 
Hamiltonians H s i(2) which are linear in sl(2) generators Y a (but with modified constants a, g) 
and have on each fixed subspace £([/«]) equidistant energy spectra obtained from Eq. (2.8a). 
However, this (quasi)equidistant approximation, in fact, corresponding to a substitution of cer- 
tain effective coupling constants g instead of true operator entities g(Y ) in Eq. (2.4), does not 
enable to display many peculiarities of models (1.1) related to essentially non-equidistant parts 
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of their spectra. Therefore, it is needed in corrections, e.g., with the help of iterative schemes 
[8,14,15]; specifically, one may develop perturbative schemes by using expansions of operator 
entities g(Y ) in Taylor series in Y as it was made implicitly for the Dicke model in [8,9]. But 
there exist a more effective, incorporating many peculiarities of models (1.1), way to amend 
the quasi-equidistant approximation. 



3 SL(2) energy functionals and variational schemes for 
solving spectral and evolution tasks 

This way is in applying SL(2) GCS |[£j];u;£) = SV(£) I [h]] v j ) as trial functions in the vari- 
ational schemes of determing energy spectra and quasiclassical dynamics [5,15]. Indeed, the 
results (2.8) are obtained by using a variational scheme determined by the stationarity condi- 
tions 

a) g"(M^O =0 , 6 ) W(W^O =0 p.!) 
ad or 

for the energy functional H([k];v;0 = ([Zi];u;£|#|[Ji];u;0 = ([h]; v\aY (Z) + Y + (£) + Y_ + 
C|[/j];f). At same time an appropriate quasiclassical dynamics, which is isomorphic to the 
exact quantum one when sl p d(2) = sl(2) [14-16], is described by the classical Hamiltonian 
equations [5,14,16] 

9 = P = ~^> n=(z{t)\H\z(t)) (3.2a) 

for "motion" of the canonical parameters p,q of the SL{2) GCS \z(t)) = exp(— z(t)Y + + 
z(t)*Y-)\ip ) as trial functions in the time- dependent Hartree-Fock variational scheme with 
the Lagrangian C = (z(t)\(id/dt — H)\z(t)); p = jcos9,q — <j>,z — 9/2exp(—i(f)) for su(2) and 
p = j cosh 0,q — (f>, z — 6/2exp(—i(f)) for su(l,l). An equivalent formulation in Y = (Y]_, Y 2 , Y ) 
space can be given in terms of sl(2) Euler-Lagrange equations, 

y = x VC, 

C = ±y 2 +y\ + yl y t = (z(t)|F^(t)), y ± = Vl ± y 2 , V = (d/d Vl ,d/dy 2 , d/dy ) (3.26) 

reducing to the well-known (linear) Bloch equations [5,14,17]. 

Similarly, general ideas of the analysis above and calculation schemes (3.1), (3.2) may be 
extended to the case of arbitrary polynomial algebras sl p d{2) by using the energy functional 
W([Zi];u;f) = ([/;]; v\H u (t)\[k\\ v) with H [k] (0 being given by Eq. (2.9). Naturally, results 
obtained in such a manner are not expected to coincide with exact solutions on all subspaces 
L([li\) due to an essential nonlinearity of Hamiltonians (2.4) and their non-equivalence (unlike 
Eq. (2.7b)) to diagonal parts of Eq. (2.9); however, they yield, evidently, most close to exact 
"smooth" (analytical) solutions (cf. [5,14]). Without dwelling on a discussion of all aspects 
of such an extension we consider in detail an application of the procedure (3.1) to the most 
wide-spread class [11] of Hamiltonians (2.4) with the su(2) dynamic symmetry which includes 
the model (2.3). 

Note that the condition (3.1a) gives e t6 = g/\g | as in the linear case, and, due to the form of 
trial functions, it is sufficiently to solve Eq. (3.1b) only for finding ground states \ v — 0;£). 
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Then, expanding r.h.s. of Eq. (2.5) in Y a power series and taking into account defining relations 
for the su(2) algebra one gets after some algebra the following expressions 

E S V U(2 \U, 6) = H([li];v; 6) = C + a(l + j) + a(-j + v) cos 2r - 2\g\ ]T E+(r; j; v), 

V) " ^ (r ' J ' ° K 2 Sm2r) (2j>!(/ - «)!(/ + 1 - v)l* 

F(-v, -v + 2j + 1; / - v + 1; sin 2 r)F(-u, -u + 2j + 1; / - u + 2; sin 2 r), 

-m r 47 \ (tanr) 2 /+ 1 (2j)! u - , ■ JA y>(Z + 1 + /) 

£j(r; j; 0) = (cos r ) (/ )i( 2j - _ / _ 1)1 V + /)> 0n- 2 (-J + /) = (2j _ /)(/ + 1) > 

(3.3) 

with F(...) being the Gauss hypergeometric function [18], for energy eig envalues^ 2 )([^];Co = 
rg/\g\) where diagonalizing values of the parameter r are determined from solving the algebraic 
equation 

= E {2j m Q ~ ^ ~ + + ^IV^H + /)}, « = -tanr (3.4) 
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For the case of the su(l, 1) dynamic symmetry Eqs. (3.3), (3.4), retaining their general structure 
form, are slightly modified due to differences in the definition (2.5) of «Sy(£) for su{2) and 
su(l, 1). Let us make some remarks concerning this result. 

1) As is seen from Eq. (3.3), its general structure coincides with the energy formula given by 
the algebraic Bethe ansatz [6], and spectral functions Ej(r;j;v) are non-linear in the discrete 
variable v labeling energy levels that provides a non-equdistant character of energy spectra 
within fixed subspaces £([/«]) at <i([Zj]) > 3. Besides, due to the square roots in expressions for 
these functions different eigenfrequencies to^ u ^ = E^ u ^ /% are incommensurable: mw™' 2 ' ^ 
nu^ 2 ^ that is an indicator of an origin of collapses and revivals of the Rabi oscillations [2,8] as 
well as of pre-chaotic dynamics [19]. Note that this dependence is impossible to get by using 
GCS related to uncoupled subsystems. 

2) The r.h.s. of Eq. (3.4) is a polynomial of the degree 2j + 1 = rf([/j]), and, in general, Eq. 
(3.4) may have 2j + 1 different roots r\ corresponding to 2j + 1 different stationary values of 
the energy functional H([li\; v ; £). Therefore, one may assume that it is possible to get more 
simple expressions for Ef(r;j; v) with any v using Ef(r;j; 0) with different roots rj. Note that 
this conjecture is valid for little dimensions d([li\) when Eqs. (3.3)-(3.4) give exact results. 
Another way to modify and to simplify the results above is in using different properties, in- 
cluding integral representations, of the hypergeometric functions F(a, b; c; x); specifically, using 
relations between hypergeometric functions [18], one can express spectral functions Ef(r;j;v) 
in terms of the hypergeometric functions 4-F 3 (...; 1) (which are proportional to the sl(2) Racah 
coefficients). 

3) Evidently, Eq. (3.3) generalizes Eq. (2.8a) for the (quasi) equidistant approximation 
abovementioned. Indeed, when replacing functions <f) n -2{~j + /) by their certain "average" 
values, series in (3.3), (3.4) are summed up, and Eq. (3.3) is reduced to Eq. (2.8a); Taylor 



series expansions of functions A/0 n _ 2 (— j + /) provide perturbative corrections related to higher 
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degrees of the an-harmonicity of Hamiltonians (2.4). Furthermore, we can get an intermedi- 
ate approximation for energy spectra if replacing in Eqs. (3.1) the exact energy functional 
7~L([li];v;£) = ([k]; v\H[ ti ] v) by their mean-field (corresponding to the Ehrenfest theorem) 
approximations 

n mfa ([k};v;0 = a < F (0 > + < Y + (0 > g(< *o(0 >)+£+(< *o(0 >) < ^-(0 > +C, 

<Y a (£)>={[%v\Y a (Z)\[lM (3.5) 

Then Eqs. (3.3)-(3.4) are very simplified retaining their main characteristic features. For 
example, for the model (2.3) we find 

C + a(l +j) + a(-j + v) cos2r - 2\g\(j - v) sin2ry / (-j + v) cos2r + j + \k\ + 1 (3.6a) 
where r is determined from the equation 

o . „ n fZ . . 9 m 7 ?'sin 2 2r .„ 

— -sm2r = cos 2ry 2j sm 2 r + \k\ + H (3.66) 

2 I#I ' 2 v /2jsin 2 r+ |ife| + 1 

(Similar expressions can be found for the point-like Dicke and the second harmonic generation 
models.) Besides, substituting Eq. (3.5) in Eqs. (3.2) one may get a mean-field approximation 
for dynamics equations reducing in the Y space representation to non-linear Bloch equations 
(cf.[5,ll]) obtained from Eqs. (3.2b) by the substitution 

xju = ([ g + g*][Myo)] 1/2 , [g - g*][Myo)] 1/2 , a + l -[ g { yi + y 2 ) + g *( yi - y 2 )]lMyoT 1/2 ), 

VC = 2(y 1 ,y 2 ,y ), Mvo) = Vo + J + \k\ + 1 (3.7) 
4) Finally, Eqs. (3.3) and (3.6a) can be used for obtaining appropriate approximations 

U s H u(2)/mfa (t) = £ Sy^oV exp( *- ) |[^><[^| Sy(^) (3.8) 

[k],v 

for the evolution operators which are transformed to the form (1.4a) with the help of the 
standard group-theoretical technique [20]. 

4 Conclusion 

So, we have obtained new approximations for energy spectra and evolution equations of mod- 
els (1.1) by means of using the mapping (2.1) and the variational schemes (3.1), (3.2) with 
the SL(2) GCS as trial functions. They may be called as a "smooth" sl(2) quasiclassical ap- 
proximations since they, in fact, correspond to picking out "smooth" (analytical) s/(2) factors 
exp(^o^+ — QY_) in exact diagonalizing operators >S(£) and in the evolution operator £/#(£). 
These approximations may be used for calculations of evolution of different quantum statistical 
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quantities (cf. [8,14]) and for determining bifurcation sets of non-linear Hamiltonian flows in 
parameter space (cf. [5]). 

Further investigations may be related to a search of suitable multi-parametric specifications 
of exact diagonalizing operators S(£) = S([£o, £1, £ 2 , •••]) using exp(£o^ + — as initial ones 

in iterative schemes which are similiar to those developed to examine non-linear problems of 
classical mechanics and optics [21] or as "principal" factors in the diagonalization schemes 
like (2.7) for Hamiltonians (1.1). From the practical point of view an important question is 
to get estimations of accuracy of approximations obtained and to make comparisons of their 
efficiency with other approximations (e.g., given in [8,9,11]). For the model (2.3) (and other 
ones with the structure polynomial ip3(x) of the third degree) it is of interest to compare results 
of approximations found above with exact calculations obtained by considering solvable cases of 
models under study. One of latters is given by integral solutions [12] and other may be yielded 
by the Riccati equations arising from a differential realization of sl p ^{2) generators V a [13]: 

\/_ = d/dz, V = zd/dz + l , V + = ip n (zd/dz + loXd/dz)- 1 (4.1) 

which is, in turn, related to a realization of sl p d(2) generators V a by quadratic forms in sl(2) 
generators Y a (cf. [15,22]). (In fact, this realization was used implicitly for obtaining exact 
integral solutioms [12].) Besides, it is also of interest to investigate possible connections of these 
results with quasi-exactly solvable sl(2) models [23,24]. 
The work along these lines is now in progress. 
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